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Abstract 

The  equation  of  motion  for  the  polarization  vector  for  a 
double  heterostructure  of  polar  cry. "als  is  solved  exactly  within  the 
framework  of  the  continuum  model.  There  exist  only  two  types  of  phonon 
modes,  the  interface  modes  and  the  confined  bulk  modes,  whose  eigenvectors 
are  obtained  explicitly.  Dispersion  relations  are  derived  analytically  for 
the  interface  modes,  while  the  confined  bulk  modes  are  dispersionless,  a 
fact  consistent  with  the  model.  It  is  also  found  that  in  Raman  scattering 
experiments,  the  symmetric  interface  modes  are  predominantly  longitudinal 
(L0)  and  the  antisymmetric  interface  modes  transversal  (TO) .  In  the  central 
region  of  the  Brillouin  zone,  however,  they  both  split  into  two  branches 
oscillating  at  L0  and  TO  frequencies,  respectively.  Possible 
reinterpretation  of  various  experiments  are  briefly  discussed. 
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I .  Introduction 

There  has  been  great  interest  in  recent  vears  in  the  studv  of  various 
vibrational  modes  supported  bv  semiconductor  heterostructures.  The  patterns 
of  normal  modes  of  vibrations  that  determine  the  electronic  properties  in 
such  structures  are  evidently  different  from  those  in  the  bulk.  The 
presence  of  interfaces  necessarily  alters  the  phonon  modes,  and  possibly 
even  their  interaction  with  electrons  may  be  modified  because  of  the  reduced 
dimensionality . 

Experimental  investigations  of  the  spectra  of  possible  acoustic  modes 
1-4 

in  superlattices  have  shown  fairly-good  agreement  with  theory.  On  the 

other  hand,  the  longitudinal  optical  (LO)  modes  in  polar  crystals  are  much 

less  well  understood  in  these  contexts.  In  all  the  calculations  such  as 

5-7  8 

free-carrier  absorption  of  light,  scattering  rates  and  polaronic 
9 

effects,  etc.,  the  usual  bulk  Frolich  Hamiltonian  is  assumed  for  the 

electron-phonon  interaction  in  confined  systems.  The  only  requirements  are 

that  the  material  elastic  properties  match  at  the  interfaces  and  that  the 

dielectric  properties  are  equal.  Mo’  locently,  the  surface  optical  (SO) 

mode  has  been  included  in  the  treatment  or  ,  olarons  confined  in  a  slab^  or 

11  12 

near  the  interface  in  semi  -  inf  inite  systems.  ’  For  these  surface  and 
interface  situations,  the  bulk  Frolich  Hamiltonian  is  still  employed  for 
treating  the  LO-phonon  contribution. 

On  the  other  hand,  evidences  of  confined  modes  peculiar  to  different 
types  of  layered  structures  nave  been  noticed  in  various  experiments. 

Measurements  of  magneto-absorption  and  cyclotron  resonance  in  GalnAs/InP 

13  14 

supe r  1  a 1 1 ic e s  and  GaAs/GaAlAs  heterostructures  indicate  that  the 


electron-LO-phonon  interaction  in  these  structures  can  be  fundamentally 
different  from  that  in  the  bulk  case.  In  a  numerical  study  of  possible 


modes  of  Che  optical  phonon  in  layered  polar  crystals,  it  is  found  that 

phonon  modes  tend  to  be  confined  in  each  layer  and  that  the  penetration  of 

vibrations  into  the  adjacent  layer  is  negl igib le . ^ ^ ^  Moreover,  the 

existence  of  confined  phonon  modes  has  been  directly  observed  in  a  GaAs 

single  quantum  well  of  GaAs/GaAlAs  heterostructures .  ^ 

In  theoretical  investigations  of  the  vibrational  modes  in  an  ionic 
1 8 

slab,  Fuchs  and  Kliewer  have  found  the  bulk  LO  mode  with  the  nodes  at  the 

surfaces  as  well  as  the  SO  modes  of  different  symmetries.  Interface  modes 

19 

have  been  derived  by  Wendler  by  considering  the  polarization  field  in  a 

20 

double  layer  system,  and  by  Lassnig  using  the  energy-loss  method  in  a 

double  heterostructure  (DHS)  of  polar  sericonductors .  An  alternative 
2 1 

treatment  of  the  DHS  predicts  some  peculiar  phonon  modes  that  have  not 
been  borne  out  by  observation. 

In  this  article  we  present  solutions  for  optical  phonon  modes  in  a 

22 

semiconductor  DHS  using  the  continuum  model  of  Born  and  Huang.  The  method 

18  19  23 

of  solution  has  been  developed  by  various  authors.  ’  '  Apart  from  the 

long-wavelength  limit  in  the  model,  no  further  approximation  is  made 
throughout  our  calculation.  Dispersion  relations  and  eigenvectors  for  all 
the  normal  modes  of  lattice  vibration  are  derived  analytically.  It  is  found 
that  there  exist  two  types  of  phonon  modes,  the  interface  modes  and  confined 
bulk  modes . 

While  the  existence  of  interface  modes  has  been  well  recognized 

experimentally ,  ^ ^  their  eigenvectors  and  dispersion  relations  in  a  DHS  P°r_ 

I 

are  solved  explicitly  for  the  first  time  in  this  paper.  Our  results  show 
that  either  the  symmetric  or  the  antisymmetric  interface  modes  has  two  '  -  n- 
branches.  Their  frequencies  at  the  center  of  the  Brillouin  zone  are  exactly  _ _ 


the  same  as  those  of  the  bulk  LO  and  TO  phonon  in  each  material. 


i  ..i.il/ii  itv  Codes 
and/or 

Dial  Special 


□  □ 


Expe  r  mie n t a  1  evidence  of  confined  bulk  modes  has  recently  been 
17  25  26 

reported.  J  '  We  find  that  both  the  bulk  LO  and  transverse  optical  (TO) 
modes  are  strictly  confined.  Further  investigation  on  implications  of  such 
confinements  is  being  carried  out  and  will  be  reported  elsewhere. 

In  Sec.  II,  we  outline  the  procedure  for  deriving  the  equation  of 
motion  for  the  polarization  vector.  The  coupled  integral  equations  are 
solved  for  the  interface  modes  in  Sec.  Ill  and  the  confined  bulk  modes  in 
Sec.  IV.  Consequences  and  implications  of  our  results  are  discussed  in  Sec. 
V. 

1 1  ■  Equation  of  motion  of  the  polarization  field 

Consider  a  DHS  of  two  different  polar  crystals  as  shown  in  Fig.  1.  A 

layer  of  material  1  with  thickness  a  is  sandwiched  between  two  thick  layers 

of  material  2.  We  take  the  z-axis  to  be  perpendicular  to  the  interfaces 

which  are  located  at  z  -  0  and  z  -  a,  respectively.  Following  Born  and 
22 

Huang,  we  start  in  the  continuum  approximation  with  the  equation  of  motion 
for  the  relative  displacement  u(r,t)  of  the  ion  pair  in  material  v  (i/  - 

1.2), 


VVr’c) 


2  -» 
-LIU)  U 
1/01/1/ 


(r,t) 


+  e*E(r,t) 

i/ 


(1) 


where  u  is  the  reduced  mass  of  the  pair  of  ions,  luo  is  the  short-range 

o 

— ►  -*♦  m 

force  constant  not  including  Coulomb  fields,  E(r,t)  is  the  local  electric 
field,  and  e*  is  the  effective  charge  of  the  ions.  The  subscript  v  labels 
the  material  considered.  The  oscillating  ions  produce  a  polarization  field 


P(r,t)  given  by 


(2) 


P(r,t)  =  n  e*u(r,t)  +  n  a  E(r,t) 

1/1/  1/1/ 

where  n  is  the  numbei  of  ion  pairs  per  unit  cell  and  a  is  the 

polarizability.  The  first  term  in  (2)  represents  the  contribution  of  the 

ion  pair  when  the  lattice  vibrates,  and  the  second  term  is  the  electronic 

polarization  of  the  ions  due  to  the  electric  field  associated  with  the 

optical  modes.  The  part  of  the  polarization  produced  by  the  electron  itself 

as  it  moves  through  the  crystal  is,  however,  not  included  in  our 

23 

consideration,  since  the  continuum  model  is  not  valid  for  such  an  effect. 
The  local  field  in  (2)  is  related,  in  the  long-wavelength  limit,  to  the 
polarization  by 

dr'  r(r-r')*P(r')  (3a) 


E(r,t)  - 


E/r.t) 


+  4ir 


I 


where 


4x 


E^(r , t)  =  —  P(r , t) 


(3b) 


and  F  denotes  the  Green  tensor  with  components 


i_ 


af}  4jt  3x  3xfl  i-» 

a  0  |  r  -  r 


(3c) 


The  equation  of  motion  for  the  polarization  then  follows  by  plugging 


Eqs .  (2)  and  (3)  into  (1)  : 


0 


;  1  -  ~a  n  ]  P  (  r  ,  t )  + 
3  v  v 


o  n  e* 

[w„..  -  4ir(a  n  t-  ^  u  ~  )  ]  P(r  ,  c) 


01/ 


1/  1/  oi/  3q 


=  47TQ  n 
1/  1/ 


dr  r(r-r' ) *P(r' ,d) 


2  n  e* 

+  4*(<W‘>o»/  +  )  |  dr'  r ( r- r '  )  •  P ( r '  ,  t )  .  (4) 


I 


The  time  -  dependent  part  of  the  polarization  can  be  separated  bv  assuming 
-*■  -*  -*  -*  Icot 

P(r,t)  -  P(r)e  which,  after  substituting  into  (4),  yields  the  equation 
for  P(r) , 


A  -X 

_ v  01/ _ 

an  (A  -A  )-l 
u  /  i/  01/ 


4tt  , 


I 


P(r)  -  4 tc  dr'  T(r-r ' )  *P(r'  ) 


(5a) 


where  we  have  defined  the  parameters 


A  -  4ttw  /w  (5b) 

v  pi/ 


01/ 


2  2 

4»(i’  /u> 

oi/  pi/ 


(5c) 


7  2 

with  the  ion  plasma  frequency  u>  -  4?m  e*  /a  . 

pi/  i/i/  i/ 

Since  the  translational  invariance  in  the  z -direction  is  destroyed  by 

the  interfaces,  we  introduce  the  two-dimensional  vectors  k  and  p  so  that  k  - 
— *  -+ 

(*,q)  and  r  -  (p,  z)  .  The  two-dimensional  Fourier  transforms  can  now  be 


written  as 


P(r) 


t>2 


elK' P  P(Z.z) 


(6) 


r  “  2)r  J  /c  expU^.p-Kjzj ) 


(7) 


Differentiating  (7)  twice  with  respect  to  the  coordinates,  we  obtain 


r  -  -  — 


l 

4* 


■ 


cbc  e 


\K*  f>  -  K  Z 


KK/2  7r« 


(8) 


where 


K  -  (>c ,  i9(z)*) 


(9) 


with  the  step  function  9(z)  -  ±1  for  z  <  0.  Substituting  (6)  and  (8)  into 
(5a),  and  moving  the  term  -P  from  the  right-hand  side  to  the  left,  we  can 
wrice  afcei  ^>ou,c  a igebia  che  resulting  equation  in  a  mOie  symmetric  form  as 


4* 


x*1^)  0  0 

0  X^L(W)  u 

0  0 


hit,. 


2w 

K 


L 


dz '  e 


-k(z-z' ) 


KK*  P(k , z ' ) 


(10) 


where 


*„<w) 


is  defined  by 


(w) 


A  -A 
t'  ot/ 


a  n  (A  -A  )  -  1 

1/  1/  1/  Oi/ 


4* 

3 


(11) 
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It  turns  out  that  y^  ( w )  is  the  isotropic  dielectric  susceptibility  and  is 
related  to  the  dielectric  function  bv  y  (u)  -  <r  <■_•)  -  1  with 

v  1/ 


£  (W) 


) 


^Lt/ 

2 


2 

2 


i  12a) 


e  -  1  +  4  war  n  /(I  -  ~  a  n  ) 
&!/  1/1/  1/1/ 


(12b) 


where  we  have  defined  the  LO  and  TO  phonon  frequencies 


2  2  2  2  ...  gj 

o>,  -oj  +Tw>  /(I  +  ~a  n  ) 
Li/  oi/  3  pi/  inn 


2  2  1  2...  4w  . 

OJ  -  o)  -  u  /(l  --an) 

Ti/  oi/  3  pi/  3  i/  u 


(13a) 


(13b) 


Since  the  interface  phonons  propagate  in  the  xy-plane,  it  is  more 

convenient  to  express  the  polarization  vector  as  P  -  ( w ,  P  )  where  x  is  a 

— ♦ 

two-dimensional  vector  defined  by  n  -  (P  ,P  ).  Thus 

K  Z 


P(K.z)  -  P  (k,Z)k  +  P  (*,z)z  +  P  (k,z)s  ,  (14) 

K  Z  S 

where  the  unit  vector  s  is  defined  by  s  -  z  *  k.  Substituting  (14)  in  (10), 
we  can  separate  the  s -component  and  decouple  (10)  into  two  equations: 


X  L(<w) 


0 


0 


X  *  (<*>)«  „(<->) 


*(*, z) 


tr. 


dz'M(z-z' )  •  ir(<c , z '  )  (15a) 


tor  the  s  -called  p - po lar ization ,  and 


V  1  -w)  P  !  *  ,  z  )  =  0  (  15b ) 

s  s 


for  the  s -polarization,  where  M  is  a  Hermitian  matrix  given  bv 


M(z-z')  -  M'(z'-z)  =  -2nKe 


-<c|  z-z'  I 


i8(z- 


.  (16) 


i9 ( z  -  z  '  ) 


As  we  shall  show  in  the  following  sections,  Eq  .  (15a)  defines  an 
eigenvalue  problem  whose  solutions  describe  the  interface  modes  and  the 
confined  modes  Equation  (15b)  describes  the  s - polarization,  which  is  not 
of  concern  in  the  present  paper.  The  eigenvectors  it(k.z)  form  a  complete 
orthonormal  set.  Here  we  just  give  without  proof  the  orthonormality 
relation 


ret  (  k.  ,  z )  n  .  (  k  ,  z  ) 


5  .  . 


(17) 


H 

with  rj  ‘  <  u> .  ) 

V  L 


1  / f  1  +  a  n  ( A  -A  )i.  The  completeness  relation  is  given  by 

V  V  Ol/  i/ 


w*  ( ,  z  )  n  .  (k.  ,  z  '  ) 


2 

CJ 


I  5  ( z  -  z  '  ) 


(18) 


for  the  so-called  p - polarization ,  and 


(w)P  (k,z) 

s 


0 


(15b) 


for  the  s -polarization ,  where  M  is  a  Hermitian  matrix  given  by 


M(z-z')  -  M'(z'-z)  =  -2nK.e 


-K\Z-Z' 


i9(z -z ' ) 


l_  i9(z  - z  '  ) 


.  (16) 


As  we  shall  show  in  the  following  sections,  Eq .  (15a)  defines  an 
eigenvalue  problem  whose  solutions  describe  the  interface  modes  and  the 
confined  modes.  Equation  (15b)  describes  the  s -polarization ,  which  is  not 

— f 

of  concern  in  the  present  paper.  The  eigenvectors  tt(k,z)  form  a  complete 
orthonormal  set.  Here  we  just  give  without  proof  the  orthonormality 
relation 


r 


4  4 

T)  (w.)*7  (w.) 

u  i  if  i  -*  -*  -*■  -* 

dz  ~  k*(k,z)  tt^(k,z) 

pu 


5u  • 


(17) 


with  t)  (w.)  -  1/ [  1  +  a  n  (A  -A  )].  The  completeness  relation  is  given  by 

l/l  1/1/  01/1/ 


2 

01 

—Elf 


I  5 (z - z ' ) 


(18) 
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where  I  stands  for  the  unit  matrix.  Finally,  we  note  that  the  polarization 
vector  must  be  real,  and  consequently 


P(k , z)  =  P*(-k ,z) 


(19) 


Similarly,  we  have 


r-i  o  i 


0  1 


ir*(-ic,z) 


(19b) 


HI-  Interface  Modes 

It  is  easier  to  solve  the  coupled  integral  equations  (15a)  by  first 
transforming  them  into  differential  equations  .  This  can  be  done  by 
differentiating  (15a)  with  respect  to  z  twice  and  at  the  same  time  requiring 


det 


(w)  0 


-1. 


0 


*  0 


(20) 


The  resulting  equations  are 


d_ 

dz 


ix  P  (k , z ) 
z 


(21a) 


2-*  -*■ 

-  x(k,z)  -  k  *(*,z) 
dz 


(21b) 
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The  solutions  to  (21)  take  the  form 


P  (k,z)  = 
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■  A  KZ 

iA2e 

. , .  kz  -KZ, 
iCA^e  -B^e 


-  iB^e 


•kz 


z  <  0 
0  <  z  <  a 
z  >  a 


(22a) 


,  .  KZ 

A2e 

*  Z 

< 

0 

P  (k , z)  -  - 
z 

KZ  „  -KZ 

A^e  +B^e 

,  0 

< 

z  <  a 

(22b) 

n 

1  B2e 

,  z 

> 

a 

Substituting  (22)  in  the  integral  equation  (15a),  one  obtains  a  set  of 
homogeneous  equations  for  the  amplitudes  A^  and  B^  of  the  p  -  polar  ization . 
The  condition  for  the  existence  of  a  nontrivial  solution  then  leads  to  the 
dispersion  relation 


<l(^)  -<2^> 

€1(ca))+«2(w) 


-  ±e 


ica 


(23) 


The  +  and  -  signs  on  the  right-hand  side  of  (23)  correspond  to  the  symmetric 
and  antisymmetric  modes  of  the  interface  phonons,  respectively.  The 
polarization  amplitudes  are  found  to  satisfy  the  following  relations: 


B-/A-  -  B-/A- 


+ 


xa 


e 


+ 


*2(«) 

x1(«) 


(1  ±  e*3) A^ 


(24a) 


(24b) 
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+  x2(») 


/cav 


B7  -  ^7—  (lie")  B7 
2  X, (w)  1 


(24c) 


Equations  (24)  are  equivalent  to  the  boundary  conditions  that  the  wave 
functions  have  to  satisfy  at  the  interfaces.  Thus  the  differential 
equations  (21)  yield  naturally  the  correct  boundary  conditions  after 
substituting  their  solutions  (22)  into  the  integral  equation  (15a).  On  the 
contrary,  the  hydrodynamic  terms  introduced  in  Ref.  23  are  inconsistent  with 
the  boundary  conditions.  Combining  Eqs .  (22)-(24),  we  find  the  eigenvectors 
for  the  antisymmetric  interface  phonon  modes  to  be 


jr  - 


f „ (w) - 1 

f  VT^T'  sinh<?> 

C  [i  sinhlK(z-f) )  ,  cosh{ie(z-a/2) ) 
a  z 

k( z-a) 


l  C  [ 


« 2(w) *1 


al c ^(w) -1 


]  e 


/caN 


(1,-1)  sinh(~) 


z  <  0 

0  <  z  <  a 

z  >  a 


(25a) 


and  for  the  symmetric  modes  to  be 


«2(w) 


C  [i  cosh{«(z-a/2) } ,  sinh{«(z-a/2) ) ] 
e0(w)-l  .  . 

.  _  .  2  .  -*(z-a)  . .  ..  ,  ,k  a. 

I  e  O-.-D-  cosh(2  ) 


«a. 


z  <  0 


sl€l(w)-lJ 


0  <  z  <  a 


z  >  a 


(25b) 


where,  according  to  (17),  the  normalization  constants  are  given  by 


a ,  s 


J  ~ 


k  ^1  ^2  €1  *2  2 

- - /[ — -  -  — -  —(—')  } 

sinh(*a)/l  2  2 

W  1  U  „  Z  i 

pi  p2 


(25c) 
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A  heterostructur e  composed  of  two  media  with  dielectric  functions 
given  by  (12)  always  supports  four  distinct  interface  modes  of  vibrations, 
two  from  each  medium.  The  dispersion  relations  for  these  modes  can  then  be 
calculated  explicitly  from  (12)  and  (23).  The  results  are 


{'-2(w 


22  2  2.  *a. 

T1+WL2 }  +  €®1^wT2+<<,L1  cothT) 


.,22  2  2  2  2  2.2  2,*a. 

-  1  €oo2(uTl'WL2)  £®1  WT2  "WL1^  COth  2  ^ 


+  2£®l£®2rL  (wTl+tJL2)  (‘JT2'HJLl) 


'  2(CJT2WL1+WL2‘JT1)1  COCh(?))li}li/l2t£co2  +  e®l  coth(^)])4  (26a) 

A  ■  {t-2(a,Tl+wL2)  +  £«1(WT2+WU)  tanh(f } 

.,2.2  2  2  2  2  2  ,2  ,2  .xa 

±  U-2(<JTr°JL2)  +  £®l(wT2"a,Ll  tanh  (2  } 


+  2£col£«2[(wTl+UJL)(wT2hJLl) 


2(WT2WL1+^T1“L2)1  tanh(r))4}4/(2[£»2+£«l  tanh(?)1,! 


(26b) 


It  is  seen  from  (26)  that  the  interface  phonon  energies  depend  explicitly 
on  the  dimensionless  quantity  *a. 

K,  3. 

Let  us  now  look  at  the  limiting  cases.  When  a  -»  «>,  tanh(~  )  -  1  and 

/C  3i 

coth(~)  -  1.  Therefore,  both  (26a)  and  (26b)  approach  the  same  limit, 


given  by 


± 

<jj 

a o 


) 


+  f®2  ^T1+U>L2) 


+ 


2  2  2  2 


2.2  2  2 

£®2  WL2 "^Tl 


+  2 


£=ol£»2[Kl+a,L2)(wT2+wU) 


2(“T2WL1^1WL2)]) 


x 


[2(£®i+ea.2)] 


(27) 


which  is  identical  to  the  result  of  a  bilayer  system  with  only  one 
19 

interface,  as  it  should.  In  the  limit  a  -*  0 ,  the  system  reduces  to  a  bulk 
material  2  with  frequencies  anc*  yhen  k  -*  0 ,  tanh(~)  -  0  and 

/C  <1 

co th(— )  -*■  We  then  find  from  (26)  that 


/  r  2  2  ,  ,  2  2  N  , 

y[WT2+WLl  1  (WT2-WL1)1/2 


a,T2’WLl 


(28a) 


Jl4^2L2  ±  <»n-“L2>l''2 


WT1’WL2 


(28b) 


that  is,  the  limiting  frequencies  are  given  by  the  bulk  LO  and  TO 
frequencies  of  the  two  materials.  It  may  be  worth  mentioning  yet  another 
limit  at  this  point.  When  the  characteristic  parameters  of  the  two 
dielectrics  approach  each  other,  or  when  -*  e^,  we  find  for  a  given  width 
that  the  amplitudes  of  the  interface  modes  of  vibration  diminish 
continuously  and  become  zero  at  -  e^. 

The  interface  phonon  modes  have  no  connection  with  bulk  polarization 
charges  because  V«P  -  0.  They  are  accompanied  by  the  surface  charges  a 


at  the  interfaces.  These  charge  densities  can  easily  be  determined  by 


15 


calculating  the  difference  of  the  polarization  eigenvectors  in  the  z- 
direction  on  both  sides  of  the  interface  concerned.  Thus  we  find  from  the 
z-components  of  t  in  (25) 

S  ,  cl 


(29a) 


for  the  symmetric  mode  and 


(29b) 


for  the  antisymmetric  mode,  where 


X1  sinh(— )  +  x2  cosh(— ) 


*l"l  X2*2 


w  , 

pi  pi 


,  ,/ca. 
coth(~ ) 


(30a) 


(30b) 


IV .  Confined  bulk  L0  and  TO  modes 

For  the  LO  modes,  -  uT  ,  and  from  (12a)  we  have  e  (m,  )  -  0  and 

Lu  u  Lu 

e  ,  (wT  , )  r*  0  for  u '  *  v .  Thus  D  -  0  in  layer  v ,  and  E  -  -  P/e  The 
v  Im  o 

differential  equations  satisfied  by  the  polarization  field  associated  with 
the  longitudinal  phonon  modes  follow  from  (15a) .  In  layer  v  we  have 
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4“  PL(«,z)  =  i/c  PL  (ac,z)  ,  (31a) 

dz  K  Z 

and  in  layer  v'  *  v 

2 

P^(ie.z)  -  «2P*>,z)  .  (31b) 

«  Z  A.  Z 


The  solution  to  (31b)  has  the  same  form  as  (22)  with  the  coefficients 
determined  by  the  boundary  conditions  for  the  field  vectors.  In  layer  v , 

-  0,  and  E  -  -  P/e  .  From  the  continuity  of  D  across  the  interface,  we  have 
D  =  0  in  layer  u'  .  Since  e  ,  *  0 ,  we  must  have  E  =  0  and  P  =  0.  Therefore 
(31b)  has  a  solution  that  is  identically  zero  everywhere  in  the  layer  u '  . 
The  boundary  conditions  in  layer  v  can  be  found  from  (15a)  for  -  0  and  \ 

-  -1.  After  some  algebra,  we  find  from  the  coupled  equations  that 


-*L  -*L 

P  (a)  -  P  (0)  -  0 


Hence,  Eq.  (31a)  is  satisfied  in  layer  1  by  the  eigenvectors 


(32) 


m 


0 


_L.  .  .  , m w  .  row  .row  . 

C  l  sm(  z),  cos  (  z) 
m  a  a/c  a 


z  <  0 
0  <  z  <  z 
z  >  a 


(33a) 


where  m  is  an  integer,  and  in  layer  2  by 


cu 


17 


-L 

n 

q 


C  [i  sin  (qz),  ^  cos  (qz) ] 


0 


C^[i  sin  q(z-a),  ^  cos  q(z-a) 
q  k. 


z  <  0 
0  <  z  <  a 
z  >  a 


(33b) 


where  q  is  real,  with  the  corresponding  eigenf requencies  and  w  , 
respectively.  These  modes  are  highly  degenerate  vibrations.  We  remark  that 
the  dispersionless  nature  of  the  confined  bulk  modes  is  consistent  with  the 
long-wavelength  limit  which  is  implied  in  the  continuum  model.  The  bulk 
optical  photons  have  constant  energy  as  k  -  0. 

Equations  (33)  clearly  show  that  the  longitudinal  optical  phonons  are 
completely  confined  by  the  interfaces.  In  the  central  layer  (labeled  1), 

the  confinement  leads  to  the  quantization  q  —  nnr/a,  where  m  -  0,  ±1,  ±2 . 

while  in  the  semi  -  inf inite  side  layers  (labeled  2),  the  wave  number  q 
remains  continuous.  The  state  vectors  are  normalized  according  to  (17)  with 
the  normalization  constants  given  by 


CL  / 


ra  Jr) 7  a  /  2  22.2 

1  Jk  +  m  n  /a 


(34a) 


„L  _L2  ,1  g 

q  "  ^2  *  Jk2  +  q2 


(34b) 


The  confined  LO  phonons  are  related  to  both  bulk  polarization  charges 
and  interface  polarization  charge^.  The  former  can  be  found  from  p  -  7*P 
and  the  latter  follows  from  tne  boundary  conditions  at  the  interfaces.  The 
results  are 


Pq  “  Cq^'C  +  sin(qz) 


jL  -  C 
m  m 


nm  f(-l) 

“  Vd” 


at 


0 

a 


(35b) 


(36a) 


at 


z  —  0 


z  -  a 


(36b) 


We  now  turn  cur  attention  to  the  TO  ohonons  for  which  w  =  w_ 

Tv 

Equation  (12a)  then  implies  that  y  ^ (<j  )  =  0  and  y  ,)  0  for  i/1  *  y. 

v  iv  v  Jv 

Hence  E  -  0  and  D  -  P  in  layer  v .  The  same  consideration  and  procedures  as 
described  above  for  the  L 0  phonons  lead  to  eigenf requencies  and  with 

corresponding  eigenvectors 


■  0 

z  <  0 

-T  r  .  mi  .  nw  ,  .  .imr  .  , 

Ci  cos — z),  sin( — z) 

ra  a/c  a  a 

0  <  z  <  a 

(37a) 

[  0  z  >  a 


where  m  is  an  integer,  and 


-T 

n 

q 


C^[i^  cos (qz ) ,  sin(qz)] 

•  0 

T  Q 

-  Cq^K  COS  <3(z*a)-  sin  q(z-a)] 


z  <  0 

0  <  z  <  a  ,  (37b) 

z  >  a 


where  q  is  a  real  number.  We  see  from  (37)  that  the  TO  phonons  are  also 
strictly  confined  by  the  presence  of  interfaces.  Once  more, 


these 


eigenvectors  are  normalized  according  to  (17),  and  the  normalization 

constants  are  given  by  the  same  expressions  as  (34)  except  for  the 

replacement  of  o>T  by  w_  The  confined  TO  modes  are,  however  not 

L  u  Ti/ 

associated  with  any  polarization  charge,  neither  bulk  nor  surface  charge. 

The  s -polarization  modes  are  given  by  the  solution  of  (15b)  with 

X  ^(w)  =  0.  This  implies  that  the  s - polar ization  modes  exist  onlv  when  the 
v 

eigenfrequencies  are  those  of  the  transverse  optical  phonons  in  either 
medium.  Since  these  modes  are  completely  decoupled  from  the  other 
vibrational  nodes,  they  are  net  involved  in  the  interaction  with  electrons 
and  hp-i.co  will  p<tt  ho  discussed  f”rther. 

V .  Discussion 

We  have  shown  that  there  exist  two  types  of  phonon  modes  in  a  double 
heterostructure  consisting  of  two  semiconducting  materials,  the  interface 
phonons  and  the  confined  bulk  phonons.  The  interface  modes  may  be  either 
symmetric  or  antisymmetric  with  respect  to  the  center  of  the  system.  They 
are  dispersive  in  nature,  and  their  frequencies  for  gi  .  en  materials  depend 
solely  upon  the  dimensionless  quantity  k.  a.  In  the  center  region  of  the 
Brillouin  zone,  these  modes  have  the  same  frequencies  as  those  of  the  bulk 
LO  and  TO  phonons  in  each  material.  For  this  reason,  we  shall  refer  to  them 
as  "LO-like"  and  "TO- like"  interface  phonons. 

Since  the  bulk  frequencies  are  determined  by  the  positions  of  the 
zeroes  and  poles  of  the  dielectric  functions  as  can  be  seen  from  (12), 
different  compositions  of  the  double  heterostructure  can  result  in  different 
frequency  combinations.  However,  only  three  distinct  combinations  as  shown 
in  Fig.  2  are  possible,  where  we  have  assumed  without  loss  of 

generality.  It  is  observed  that  as  the  width  of  the  central  layer 


:o 


increases,  the  four  interface  modes  become  two  degenerate  modes.  In  the 
limit  of  large  a,  these  modes  have  the  same  frequencies  as  those  in  a 
bilaver  he  te  ro  s  t  rue  tur  e  .  In  case  (a),  the  degenerate  modes  are  material- 
like,  while  in  the  other  two  cases  they  are  LO-like  and  TO-like. 
Experimentally,  only  case  (a)  has  been  observed  thus  far.  It  is  therefore 
interesting  to  carry  out  experiments  on  samples  with  u>  w  „  >  u  w  such 

i^X  L-iZ.  IX  L  Z. 

as  GaAs/Ga^  ^a-2q  yAs  (GaAs-type)  and  InP/AiSb. 

It  should  also  be  of  great  interest  to  note  that  the  peculiar  mode 

27 

observed  in  the  90°  Raman  scattering  experiment  may  be  understood,  as  has 

2  8 

been  pointed  out  in  a  recent  communication,  in  terms  of  the  interface 

modes  derived  in  Sec.  III.  In  other  words,  the  novel  slab  modes  reported  in 

Ref.  27  are  in  fact  the  interface  modes.  A  detailed  analysis  of  this 

experiment  will  be  published  elsewhere,  and  here  we  give  only  a  qualitative 

account.  Since  Raman  scattering  experiments  involve  only  phonons  of  very 

small  k  ,  the  dominant  component  of  the  polarization  vector  n  is  P 

a  z 

according  the  (25a).  Therefore  the  antisymmetric  interface  modes  are 
predominantly  TO  modes.  Tn  the  central  layer,  this  TO  mode  oscillates  at 
the  LO  frequency  of  GaAs ,  in  agreement  with  the  experimental  result  in  the 
right-angle  scattering  configuration.  Similarly,  (25b)  shows  that  the 
polarization  has  a  dominant  P^  component,  or  the  symmetric  interface  mode 
in  the  central  layer  is  predominantly  longitudinal  and  oscillates  at  the 
bulk  TO  frequency  of  GaAs . 

In  addition  to  the  Raman  scattering  experiments,  the  interesting 

pinning  phenomenon  has  been  reported  in  recent  measurements  of  cyclotron 

resonance.  That  the  electron  interacts  with  optical  phonons  at  the  bulk  TO 

frequency  has  been  observed  in  the  measurements  of  the  magnetopolaron 

29-31 

frequency  in  semiconductor  quantum  wells,  and  it  has  been  attributed  to 


1 


29 

the  classical  dielectric  effect.  This  is  essentially  a  polariton  effect 
rather  than  a  polaronic  one.  When  the  ls-2p  transition  energy  of  a 

hvdrogenic  impurity  atom  in  a  GaAs  quantum  well  is  measured  in  strong 

31  - 1 

magnetic  fields,  the  pinning  is  found  at  a  frequency  about  40  cm  below 

(-  20  cm  ^  below  .  To  our  knowledge,  there  is  no  theory  up  to  the 

32 

present  time  that  can  account  for  this  result.  The  existence  of  traveling 

L0  phonons  and  the  zone- folding  effect  has  been  suggested  as  a  possible 

ai 

source  of  this  phenomenon.  We  have  solved  the  interface  phonon  modes  in  a 
supe r lattice ,  and  our  preliminary  results  indicate  that  probably  the 
interface  modes  are  responsible  for  this  strange  pinring  phenomenon.  More 
careful  study  is  necessary,  however,  before  any  definite  conclusion  can  be 
made.  Work  along  this  direction  is  also  underway  and  will  be  discussed  in 
forthcoming  publications. 
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Geometrv  of  Che  double  hecerostructure . 


Dispersion  relations  of  the  interface  modes  in  the  double 
heterostructu~e  for  different  compositions:  (a)  GaAs/AiAs  for  which  « 

>  l>t2  >  u>L1  >  wT1>  (b)  GaAs/AiQ  7As  (GaAs-type)  for  which 

>  wT2  >  w  (c)  InP/AiSb  for  which  >  u>T^ . 
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